The energy analysis for the monte carlo simulations of a diffusive shock 
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ABSTRACT 

According to the shock jump conditions, the total fluid's mass, momentum, and 
energy should be conserved in the entire simulation box. We perform the dynamical 
Monte Carlo simulations with the multiple scattering law for energy analysis. The 
various energy functions of time are obtained by monitoring the total particles' mass, 
momentum, and energy in the simulation box. In conclusion, the energy analysis 
indicates that the smaller energy losses in the prescribed scattering law are, the harder 
the energy spectrum produced is. 

Subject headings: acceleration of particles — solar energetic particles (SEP) — cos- 
mic ray (CR) — shock waves 



1. Introduction 

The gradual solar energetic particles with a power-law energy spectrum are generally thought 
to be accelerated by the first-Fermi acceleration mechanism at the in terplanetary shocks (IPs) 
dAxford et all Il977l : iKrvmskvl 1 19771: lielll 1 19781 : blandford & Ostrikeij[l978.) . It is weU known 
that the diffusive shock accelerated the particles efficiently by the accelerated particles scattering 
off the instability of Alfven waves which are generated by the accelerated particles themselves 
dOoslinget. al. Ill98ll : Icane et.al. l[l99ol : Lee & Rvanll 19861 : berezhko. etiaTlboOBl: IPelletier. et.al 
20061) . The diffusive shock acceleration (DSA) is so efficient that the back-reaction of the ac- 
celerated particles on the shock dynamics cannot be neglected. So the theoretical challenge is 
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acceleration, there are largely three basic approaches: stationary Monte Carlo simulations, fully 
numerical simulations, and semi-analytic solutions. In the stationary Monte Carlo simulations, the 
full particle popul ation with a prescribed scattering law is ca lculated based on the particle-in-cell 



(PIC) techniques (lEUison et al. 



1996l ; IVladimirov et al.ll2006l) . In the fully numerical simulations. 



a time-dependent diffusion-c onvection equation for the CR transport is solved wit h coupled gas 
dynamics conservation laws (IKang & JonesI 120071 ; IZirakashvili & Aharonianll2010l) . In the semi- 
analytic approach, the stationary or q uasi- stationary diffusion-convection equ ations coupled to the 
gas dynamical equations are solved ([Blasi. et. al.ll2007l ; iMalkov. et. al.ll2000|) . Since the velocity 
distribution of suprathermal particles in the Maxwellian tail is not isotropic in the shock frame, the 
diffusion-convection equation cannot directly follow the injection from the non-diffusive thermal 
pool into the diffusive CR population. So considering both the quasi- stationary analytic models 
and the time-dependent numerical models, the injection of particles into the accele ration mech 



anism is based on an assumption of the transparen cv function for thermal leakage (IBlasi. et. al 



2005l ; lKang & Jonesll2007l ; lvainio & LaitinenlboOvl) . Thus, the dynamical Monte Carlo simulations 
based on the PIC techniques are expected to model the shock dynamics ti me-dependently and also 
can e l iminate the suspicio n arising from the assumption of the injection (|Knerr. Jokipii & Ellison 
19961 ; IWang & Yanll201lh . In plasma simulation (PIC and hybrid), there is no distinction between 
thermal and non-thermal particles, hence particle injection is intrinsica lly defined by the pres cribed 
scattering properties, and so it is not controlled with a free parameter (|Caprioli. et. al.ll2010|) . 



Actually, Wang et al. IWang & YanI (|201 1[) have extended the dynamical Monte Carlo mod- 
els with an anisotropic scattering law. Unlike the previous isotropic prescribed scattering law, the 
Gaussian scattering angular distribution is used as the complete prescribed scattering law. Accord- 
ing to the extended prescribed scattering law, we obtained a series of similar energy spectrums 
with little difference in terms of the power-law tail. However, it is not clear how such a prescribed 
scattering law can affect the particles' diffusion and the shock dynamics evolution. To probe these 
problems, we expect to diagnose the energy losses in the simulations by monitoring all of the 
behaviors of the simulated particles. 

In the time-dependent Monte Carlo models coupled with a Gaussian angular scattering law, 
the results show that the total energy spectral index and the compression ratio are both effected by 
the prescribed scattering law. Specifically, the total energy spectral index is an increasing function 
of the dispersion of the scattering angular distribution, but the subshock's en ergy spectral index i s 
a digressive function of the dispersion of the scattering angular distribution (|Wang & YanI 1201 11) . 
In the dynamical Monte Carlo simulations, one find that it is the only way for the particles to 
escape upstream via free escaped boundary (FEB). With the same size of the FEB which limited 
the maximum energy of accelerated particles, we find that different Gaussian scattering angular 
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distribution generate different maximum energy particles through the scattering process at the same 
simulation time. 

In an effort to verify the efficiency of the energy transfer from the thermal to superthermal and 
the effect of the shock dynamics evolution on the shock structures, we perform a dynamical Monte 
Carlo code on Matlab with Gaussian scattering angular distribution by monitoring the particles' 
mass, momentum and energy as functions of time. Our Gaussian scattering angular distribution 
algorithm consists of four cases involving four specific standard deviation values. This aim is to 
know if the various particle's loss functions are dependent on the prescribed scattering law and the 
various kinds of losses can directly determine the total compression ratio and final energy spectral 
index with the same timescale of the shock evolutions and the same size of FEB. 

In Section[2l the basic simulation method is introduced with respect to the Gaussian scattering 
angular distributions for monitoring the particles' mass, momentum and energy as function of time 
in each case. In Section [3l we present the shock simulation results and the energy analysis for all 
cases with four assumptions of scattering angle distributions. Section |4] includes a summary and 
the conclusions. 



2. Method 



The Monte Carlo model is a general model, although it is considerably expensive computa- 
tionally, and it is important in many applications to include the dynamical effects of nonlinear DSA 
in simulations. Since the presc ribed scattering law in Monte Carlo mod el instead of the field cal- 
culation in hybrid simulations (iGiacalonel 120041: IWinske & OmidilllOllh . we assume that particle 
scatters elastically according to a Gaussian distribution in the local plasma frame and that the mean 
free path (mfp) is proportional to the gyroradius (i.e., X rg), where rg = pc/ (qB), and its value is 
proportional to its momentum. Under the prescribed scattering law, the injection is correlated with 
those "thermal" particles which manage t o diffuse the shock front for obtaining additional energy 
gains and become superthermal particles (|Ellison et al.ll2005|) . 



However, the basic theoretical limit to the accelerated particle's energy arises from the ac- 
celerated particle's La rmor radius, which must be smaller than the dimensions of the acceleration 
region by at least Vs/c jHillaslll984h . where Vs is the shock velocity. The limitation to the maximum 
energy of accelerated particles due to the large Larmor radius would be ameliorated if the scattering 
angular distribution is varied. In these simulations, the size of the FEB is set as a finite length scale 
matched with the maximum diffusive length scale. Here, we further investigate the possibility that 
the accelerated particles scatter off the background magnetic field in the acceleration region with 
not only an isotropic scattering angular distribution, but also with an anisotropic scattering angular 
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distribution. Actually, this anisotropic distribution would probably produce an important effect on 
simulation results. With the same size of FEB, the scattering law applied isotropic distribution 
or anisotropic distribution would produce different maximum energy part icles. So an anisotropic 
scattering law in the theory of the CR-diffusion is also needed (|Bellll2004|) . 



Case A (o-i/4,|ii-0 ) 



Case B (o-i/2,|ii-0 ) 





Fig. 1. — The entire evolutional velocity profiles in four cases. The dashed line denotes the FEB 
position in each plot. The precursor is located in the area between the downstream region and the 
upstream region in each case. 

The particle-in-cell techniques are applied in these dynamical Monte Carlo simulations. The 
simulation box is divid ed up into some number of cells and the field momenturn is calculated at 
the center of each cell (|Forslundlll985l: ISpitkovskyI 120031: INishikawa. et. al.ll2008|) . The total size 
of a one-dimensional simulation box is set as X^ax and it is divided into the number of grids N„rax- 
Upstream bulk speed Uq with an initial Maxwellian thermal velocity Vt in their local frame and 
the inflow in a "pre-inflow box" (PIB) are both moving along one-dimensional simulation box. 
The parallel magnetic field Bq is along the x axis direction in the simulation box. A free escaped 
boundary (FEB) with a finite size in front of the shock position is used to decouple the escaped 
particles from the system as long as the accelerated particles beyond the position of the FEB. The 
simulation box is a dynamical mixture of three regions: upstream, precursor and downstream. The 
bulk fluid speed in upstream is U = Uq, the bulk fluid speed in downstream is t/ = 0, and the 
bulk fluid speed with a gradient of velocity in the precursor region is Uq > U > 0. Because of 
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the prescribed scattering law instead of the particle's movement in the fields, the injected particles 
from the thermalized downstream into the precursor for diffusive processes are controlled by the 
free elastic scatter mechanism. To obtain the information of the total particles in different regions at 
any time, we build a database for recording the velocities, positions, and time of the all particles, as 
well as the index and the bulk speeds of the grids. The scattering angle distributions are presented 
by Gaussian distribution function with a standard deviation o, and an average value /u involving 
four cases: (1) Case A: o = 71/4, ^ = 0. (2) Case B: a = k/2, ^ = 0. (3) Case C: o = tt, = 0. (4) 
Case D: a = oo, ;U = 0. These presented simu lations are all based on a one-dimensional simulation 
box and the specific parameters are based on IWang & YanI (|201l|) . 



3. Results & analysis 

We present the entire shock evolution with the velocity profiles of the time sequences in each 
case as shown in Figure \T\ The total velocity profiles are divided into three parts with respect to 
the shock front and the FEB locations. The upstream bulk speed Uq dynamically slows down by 
passing though the precursor region, and its value decreases to zero in the downstream region (i.e. 
Ud = 0). The precursor explicitly shows a different slope of the bulk velocity and different final 
FEB locations in different c ases. The present ve locity profiles are similar to the density profiles in 



the previous simulations by IWang & YanI (|201lt) . and the different prescribed scattering law leads 



to the different shock structures. 

The various losses of the particles and the calculated results of the shocks at the end of the 
simulation for the four cases are listed in Table [T] The initial box energy is ^o- The subshock's 
compression ratio r„,/, and the total compression ratio rtot are calculated from the fine velocity 
structures in the shock frame in each case, respectively. The total energy spectral index Tjot and 
the subshock's energy spectral index r^ub are deduced from the corresponding total compression 
ratio rtot and the subshock's compression ratio rsub in each case. The Mio^s, Pioss and Eioss are the 
mass, momentum, and energy losses which are produced by the escaped particles via the FEB, 
respectively. 

We have monitored the mass, momentum and energy of the total particles in each time step 
in each case. Figure |2] shows all the types of energy functions with respect to time. The total 
energy Etot is the energy summation over the time in the entire simulation box at any instant in 
time. The box energy Ebox presents the actual energy in the simulation box at any instant in time. 
The supplement energy Epjs is the summation of the amount of energy from the pre-inflow box 
(at the left boundary of the simulation box) which enters into the simulation box with a constant 
flux. The Efeb presents a summation of the amount of energy held in the precursor region. The 
Eout indicates a summation of the amount of energy which escapes via the FEB. Clearly, the total 
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Table 1 : The Calculated Results 



Items 


Case A 


CaseB 


CaseC 


CaseD 


Mloss 


1037 


338 


182 


127 


Ploss 


0.0352 


0.0189 


0.0123 


0.0084 


Eloss 


0.7468 


0.5861 


0.4397 


0.3014 


Efeb 


0.8393 


0.5881 


0.5310 


0.5022 


F 


1.5861 


1.1742 


0.9707 


0.8036 


Etot 


3.3534 


3.4056 


3.3574 


3.4026 


Rin 


38.25% 


25.67% 


19.98% 


14.80% 


Rloss 


22.27% 


17.21% 


13.10% 


8.86% 


^tot 


8.1642 


6.3532 


5.6753 


5.0909 


^sub 


2.0975 


3.0234 


3.1998 


3.9444 


Tfo? 


0.7094 


0.7802 


0.8208 


0.8667 


^sub 


1.8668 


1.2413 


1.1819 


1.0094 


Vsh 


-0.0419 


-0.0560 


-0.0642 


-0.0733 


Vsub 


0.0805 


0.1484 


0.1613 


0.2159 


^Lmax 


11.4115 


14.2978 


17.2347 


20.5286 


Epeak 


0.1650keV 


0.1723keV 


0.1986keV 


0.2870keV 


'-'max 


1.23MeV 


1.93MeV 


2.80MeV 


4.01MeV 



The units of mass, momentum, and energy are normal to the proton mass nip, initial momentum 
Pq, and initial energy Eq, respectively. The last two rows are shown as scaled values. 

energy in the simulation at any instant in time is not equal to the actual energy in the box at any 
instant in time in each plot. 

It can be shown that the incoming particles (upstream) decrease their energy (as viewed in the 
box frame) as they scatter in the shock precursor region. If each incoming particle loses a small 
amount of energy as it first encounters the shock, this would produce a constant linear divergence 
between the curves for E^ox and Etot- Actually, the cases in these simulations produce the non- 
linear divergence between the curves for Ej^ox and Etot, consistent with Figure |2l Such behavior is 
evident from individual particle's trajectories. Physically, all the kinds of the losses occur in the 
precursor region owing to the "back reaction" of the accelerated ions and the escaped particles via 
FEB. Various energy functions obviously show that the case applying the anisotropic scattering 
law produces a higher energy loss, while the case applying the isotropic scattering law produces 
a lower energy loss. Consequently, we consider that the prescribed scattering law dominate the 
energy losses. 
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3.1. Energy losses 

With monitoring each particle in the grid of the simulation box in any increment of the time, 
the escaped particles' mass, momentum, and energy loss functions with the time are obtained and 
shown in Figure |3l Among these energy functions, the inverse flow function is obtained from the 
injected particles from the thermalized downstream to the precursor. 

Since the FEB is limited to be in front of the shock position with the same size in each case, 
once a single accelerated particle moves backward to the shock beyond the position of the FEB, we 
exclude this particle from the total system as the loss term in the mass, momentum, and the energy 
conservation equations. According to the Rankine-Hugoniot (RH) relationships, the compression 
ratio of the nonrelativi stic shock with a large Mach number is not allowed to be larger than the 



standard value of four (|Pelletierll2001|) . Owing to the energy losses which inevitably exist in the 
simulations, the calculated results show a decreasing values of the loss of the mass, momentum, 
and the energy from the cases A, B, and C to D, respectively. Simultaneously, the inverse energy 
injected from downstream to precursor is also a decreasing values of (£',>,)a = 1.5861, {Ein)B = 
1.1742, {Ei„)c = 0.9707, and (£,„)o = 0.8036 from the cases A, B, and C to D, respectively. The 
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The Escaped Mass Losses From FEB 



The Escaped Momentum Losses From FEB 




1000 1500 
Time 



The Escaped Energy Losses From FEB 





1000 1500 
Time 



The Inverse Energy Injected from Downstream to Precursor 




Fig. 3. — The four plots denote the mass losses, momentum losses, energy losses and the inverse 
energy, respectively. The solid line, dashed line, dash-dotted line and the dotted line represent the 
cases A, B, C and D in the first three plots. In the last plot, the dashed lines marked with the 
signal of dot, plus, and cross, and the solid line represent the cases A, B, C and D, respectively. 
The units are normal to the proton mass Mp, initial total momentum Pq and initial total energy Eq, 
respectively. 



accurate energy losses are the values of {Eioss)a = 0.7468, {Eioss)b = 0.5861, {Eioss)c = 0.4397, 
and {Eioss)d = 0.3041 in each case. The inverse energy is the summation of the energy loss Eio^s 
and the net energy Ef^b in the precursor (i.e. £,„ = Ejoss + Efet)- So it is no wonder that the 
total shock ratios are all larger than four because of the existence of energy losses in all cases. 
Therefore, the difference of the energy losses and the inverse energy can directly affect all aspects 
of the simulation shocks. 



3.2. Subshock structure 

Figure m shows the subshock structure in each case at the end of the simulation time. The spe- 
cific structure in each plot consists of three main parts: precursor, subshock and downstream. The 
smooth precursor with a larger scale is between the FEB and the subshock's position Xsuh, where 
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the bulk velocity gradually decreases from the upstream bulk speed Uq to Vsubi ■ And the size of 
the precursor is almost equal to the diffusive length of the maximum energy particle accelerated by 
the diffusion process. The sharp subshock with a shorter scale only spans three-grid-lengthsb in- 
volving a deep deflection of the bulk speed abruptly decreasing from Vsubi to Vsubi^ where the scale 
of the three-grid-length is largely equal to the mean free path of the averaged thermal particles in 
the thermalized downstream. The value of the subshock's velocity can be defined by the difference 
value of the two boundaries of the subshock. 



Vsub = \Vsub\— Vsubl | • ( 1 ) 

With the upstream bulk speed slowing down from Uq to zero, the size of the downstream 
region is increasing with its constant shock velocity v^/, in each case, and the bulk speed is t/ = 
owing to the dissipation processes which characterize the downstream. The gas subshock is just an 
ordinary discontinuous clas sical shock embedded in the total shock with a comparably larger scale 



( Berezhko & Ellison 111999 ) 



According to the subtle shock structures, the shock compression ratio can be cataloged into 
two classes: one class presents the entire shock named the total compression ratio Vtot and the other 
class characterizes the subshock named the subshock's compression ratio r^ub- The values of the 
two kinds of compression ratios can be reduced by the following formulas, respectively. 



rtot = ui/u2, (2) 

rsub = (Vsub + \Vsh\)/\ Vsh I , (3) 

where ui = uq + \vsh\ , ^2 = \vsh\^ and ^1(^2) is the upstream (downstream) velocity in the shock 
frame, Vsub is the subshock's velocity determined by Equation [H and the shock velocity v^h at the 
end of the simulation is decided by the following: 

^sh ~ i-^max -^sh) /Tmaxt (4) 

where X^ax is the total length of the simulation box, T^ax is the total simulation time, and X^h is 
the position of the shock at the end of the simulation. The specific calculated results are shown 
in Table [U The values of the subshock's compression ratios are {r sub) A='^ -091 5 , (r^„i,)s=3.0234, 
(rsM&)c=3-1998 and (r„,i,)/)=3.9444 corresponding to the cases A, B, C and D, respectively. The 
total shock compression ratios with the values of (rfof)^=8.1642, (rfof)fi=6.3532, (r,o/)c=5.6753, 
and {ft ot)D=5 -0909 also correspond to the cases A, B, C and D, respectively. In comparison, 
the total shock compression ratios are all larger than the standard value four and the subshock's 
compression ratios are all lower than the standard value four. Additionally, the value of the total 
shock compression ratio decreases from Cases A, B, and C to D, while the subshock's compression 
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Fig. 4. — Final subshock fine structures in the four cases. The vertical solid and dashed lines 
indicate the positions of the shock front and subshock in each plot, respectively. The horizontal 
solid, dashed, dash-dotted and dotted lines show the values of the shock velocity Vsh, subshock 
velocity Vsubi, subshock velocity Vsubi and initial bulk velocity Uq, respectively. Three vertical 
blocks in each plot represent the three deflections of velocity: precursor region, subshock region 
and downstream region. All values of the velocity are based on the box frame. 

ratio increases from Cases A, B, and C to D. These differences are naturally attributed to the 
different fine subshock structures. 



3.3. Maximum energy 



We select some individual particles from the phase-space-time database recording the all par- 
ticles' information. The trajectories of the selected particles are shown in Figure |5l One of these 
trajectories clearly shows the route of the maximum energy of accelerated particle which undergoes 
the multiple collisions with the shock front in each case. The maximum value of the energy is ap- 
parently different in each case with increasing values of {VLmax)A = 11.4115, {VL„jax) b = 14.2978, 
{VLmax)c = 17.2347, and {VL,nax)D = 20.5286 from the Cases A, B, and C to D, respectively. And 
those particles with a value higher than the cutoff energy are unavailable owing to their escaping 
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Fig. 5. — The individual particles with their thermal velocity in the local frame vs their positions 
with respect to time in each plot. The shaded area indicates the shock front, the solid line in the 
bottom plane denotes the position of the FEB in each case, respectively. Some irregular curves 
trace the individual particle's trajectories near the shock front with time. The maximum energy of 
accelerated particles in each case is marked with the value of the velocity, respectively. 

from the FEB. The statistical data show the number of escaped particles in each case decreases with 
the number of particles {nesc)A = 1037, {nesc)B = 338, {nesc)c = 182, and {nesc)D = 127 from Cases 
A, B, and C to D, correspondingly. Except for the maximum energy of the particle in each case, 
the other particles show that some of them obtained finite energy accelerations from the multiple 
crossings with the shock and some of them do not have additional energy gains owing to their lack 
of probability for crossing back into the precursor due to their small diffusive length scale. The sta- 
tistical data also exhibit that the inverse energy injected from the downstream back to upstream is 
characterized by a decreasing reflux rate of = 38.25%, {Rinjs = 25.67%, {Rin)c = 19.98%, 

and {Rin)D = 14.80% in corresponding Cases A, B, C, and D. With the decrease of the inverse 
energy from Cases A, B, and C to D, the corresponding energy losses are also reduced at the 
rate of {Rioss)a = 22.27%, = 17.21%, (Riossk = 13.10%, and {Rioss)d = 8.86% in each 

case, respectively. Although the maximum energy of the accelerated particles should be identical 
because of the limitation of the same size of the FEB in four cases, the cutoff energy values are 
still modified by the existence of the energy losses in the different cases applied with the different 
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Fig. 6. — The two plots present the final energy spectrums on the downstream and the precursor 
region, respectively. The thick solid line with a narrow peak at E =1.3105keV in each plot rep- 
resents the same initial Maxwell energy distributions in each case. The solid, dashed, dash-dotted 
and dotted extended curves with the "power-law" tail present the energy spectral distributions 
corresponding to Cases A, B, C and D, respectively. All these energy spectrum distributions are 
plotted in the same shock frame. 

prescribed Gaussian scattering laws. 



3.4. Energy spectrum 

The energy spectrums with the "power-law" tail are calculated in the shock frame from the 
downstream region and the precursor region at the end of simulation, respectively. The same 
initial Maxwellian distribution in each case is shown in each plot in Figure |6l As shown in Figure 
[6l the calculated energy spectrums indicate that the four extended curves in the downstream region 
with an increasing value of the central energy peak {E peak) A = 0.1650fcey, {EpeakjB = 0.n23keV, 
{Epeakjc = 0.l9S6keV, and {Epeakjo = O.lSlOkeV characterize the Maxwellian distributions in 
the "heated-downstream" from Cases A ,B, and C to D, respectively. The value of the total energy 
spectral index (r,o,)A = 0.7094, (r,,,f)s = 0.7802, {rtot)c = 0.8208, and {rtot)D = 0.8667 in each 
case indicates the Maxwellian distribution in the "heated-downstream" with a decreasing deviation 
to the "power-law" distribution from Cases A, B, and C to D, correspondingly. But the value 
of the subshock's energy spectral index (r„,/,)A = 1.8668, {rsub)B = 1.2413, {rsub)c = 1-1819, 
and (rsub)D = 1.0094 present the energy spectrum distribution with a "power-law" tail in each 
case implying there is an increasing rigidity of the spectrum from the Cases A, B, and C to D, 
respectively. The cutoff energy at the "power-law" tail in the energy spectrum is given with an 
increasing value of {Emax)A='i--23 MeV, (£„mx)b=1-93 MeV, {Einax)c=^-^0 MeV and (£'„mx)z)=4.01 
MeV from the Cases A, B, C and D, respectively. In the precursor region, the final energy spectrum 
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Fig. 7. — Two plots show the correlation of the compression ratio vs the energy losses and the 
correlation of the energy spectral index vs the inverse energy, respectively. The triangles represent 
the total compression ratios and the total energy spectral index of the all cases in each plot. The 
circles indicate the subshock's compression ratios and the subshock's energy spectral index of all 
cases in each plot. 

is divided into two very different parts in each case. The part in the range from the low energy to 
the central peak shows an irregular fluctuation in each case. The irregular fluctuation indicates 
that the cold upstream fluid slows down and becomes the "thermal fluid" by the nonlinear "back 
reaction" processes. And the other part in the range beyond the central peak energy shows a smooth 
"power-law" tail in each case. 

As shown in Figure |7l the two kinds of shock compression ratios are both apparently depen- 
dent on the energy losses with respect to these four presented simulations. As viewed from Cases 
A, B, and C to D, the total compression ratio is a decreasing function of the energy losses and each 
value is larger than the standard value four, the subshock's compression ratio is an increasing func- 
tion of the energy losses and each value is lower than four. However, both the total compression 
ratio and the subshock's compression ratio approximate the standard value of four as the energy 
loss decreases. According to the DSA theory, if the energy loss is limited to be the minimum, 
the simulation models based on the computer will more closely fit the realistic physical situation. 
Additionally, the energy spectral index is also fairly dependent on the inverse energy from the 
thermalized downstream region into the precursor region. 



4. Summary and conclusions 

In summary, we performed the dynamical Monte Carlo simulations using the Gaussian scat- 
tering angular distributions based on the Matlab platform by monitoring the particle's mass, mo- 
mentum and energy at any instant in time. The specific mass, momentum and energy loss functions 
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with respect to time are presented. A series of analyses of the particle losses are obtained in the 
four cases. We successfully examine the relationship between the shock compression ratio and the 
energy losses, as well as verify the consistency of the energy spectral index with the inverse energy 
injected from the downstream to precursor region in the simulation cases which are applied with 
the prescribed Gaussian scattering angular distributions . 

In conclusion, the relationship of the shock compression ratio with the energy losses via FEB 
verify that the energy spectral index is determined by the inverse energy function with time. In 
fact, these energy losses simultaneously depend on the assumption of the prescribed scattering 
law. As expected, the maximum energy of accelerated particles is limited by the size of the FEB 
according to the maximum mean free path in each case. However, there is still a fairly large 
difference between the maximum energy of the particle from the different cases with the same 
size of the FEB. We find that the total energy spectral index increases as the standard deviation 
value of the scattering angular distribution increases, but the subshock's energy spectral index 
decreases as the standard deviation value of the scattering angular distribution increases. In these 
multiple scattering angular distribution simulations, the prescribed scattering law dominates the 
energy losses and the inverse energy. Consequently, the case of applying a prescribed law which 
leads to the minimum energy losses will produce a harder subshock's energy spectrum than those 
in the cases with larger energy losses. These relationships will drive us to find a newly prescribed 
scattering law which leads to the minimum energy losses, making the shock compression ratio 
more closely approximate the standard value of four for a nonrelativistic shock with high Mach 
number in astrophysics. 

The authors would like to thank Doctors G. Li, Hongbo Hu, Siming Liu, Xueshang Feng, and 
Gang Qin for many useful and interesting discussions concerning this work. In addition, we also 
appreciate Profs. Qijun Fu and Shujuan Wang, as well as other members of the solar radio group 
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